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DATA  REDUCTION  PROCEDURES  FOR  SEA  KING 
HELICOPTER  FLIGHT  TRIALS 


by 


N.  E.  GILBERT 


SUMMARY 


The  data  reduction  procedures  used  in  obtaining  fully 
processed  data  from  raw  flight  data  for  trials  on  a  Sea  King  Mk.50 
helicopter  are  given.  The  procedures  allow  various  corrections 
and  calibrations  to  be  applied,  removal  of  noise,  and  calculation 
of  many  additional  quantities,  some  of  which  are  used  for  kine¬ 
matic  consistency  checking  purposes.  Examples  are  provided  on 
the  running  of  the  various  computer  programs  developed.  To 
assist  in  the  use  of  the  data  for  validation  of  the  Sea  King 
mathematical  model,  output  is  obtained  in  a  form  allowing  ready 
comparison  between  trials  and  model  results  on  the  same  graphs. 


POSTAL  ADDRESS:  Chief  Superintendent,  Aeronautical  Research  Laboratories, 
p.O.  Box  4331,  Melbourne,  Victoria,  3001,  Australia. 
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1. 


INTRODUCTION 


Flight  trials  were  conducted  in  1979  by  the  Aeronautical 
Research  Laboratories  (ARL)  on  a  Royal  Australian  Navy  Sea  King  Mk.50 
helicopter  [1] .  The  primary  aim  of  the  trials  was  to  record  data 
which  could  be  used  to  validate  a  mathematical  model  developed  by 
ARL  12,  3).  The  data  acquisition  system  installed  allowed  32  channels 
of  data  to  be  recorded  on  magnetic  tape  in  serial  digital  form  at  a 
maximum  sampling  frequency  of  60  Hz  for  any  channel  [4] .  The  purpose 
of  this  document  is  to  describe  the  data  reduction  procedures  used 
in  obtaining  fully  processed  data  from  raw  flight  data. 

Section  2  describes  how  data  for  each  flight  record  are 
converted,  using  a  transcriber  unit  [5]  and  a  transcription  computer 
program,  to  a  disk  file  form  to  enable  ready  processing  on  the  DEC 
system- 10  computer  at  ARL.  The  complete  data  reduction  procedure, 
described  in  Section  3,  may  be  divided  into  several  stages.  These 
are  input  and  basic  processing,  correction  for  drop-outs,  correction 
for  instrument  and  analogue  filter  time  delays,  conversion  to 
engineering  units,  digital  filtering,  smoothing,  calculation  of 
additional  quantities,  and  output.  Some  of  the  additional  quantities 
calculated  are  for  kinematic  consistency  checking  purposes  [6] .  This 
allows  the  degree  of  confidence  in  the  flight  data  to  be  assessed  and 
suitable  corrections  to  instrument  time  lags  and  calibration  constants 
to  be  made.  Since  reference  is  made  to  the  axes  systems  in  which  both 
measured  and  calculated  quantities  are  specified,  these  axes  systems 
and  their  relative  orientation  are  first  defined  in  Fig.  1. 

Examples  of  the  running  of  the  various  processing  computer 
programs  form  the  contents  of  Section  4.  These  examples  show  how  (a) 
information  (e.g.  calibration  factors  and  offsets)  is  specified  on 
the  way  in  which  each  channel  is  processed,  (b)  the  data  reduction 
procedure  is  checked  for  correctness  step-by-step,  and  (c)  the  fully 
processed  data  is  output  in  various  forms.  Because  the  same  general 
purpose  output  program  [ 7 ,  8]  is  used  for  data  obtained  from  both  the 
flight  trials  and  mathematical  model ,  comparisons  on  the  same  graphs 
jure  readily  made. 

2 •  TRANSCRIPTION  OF  FLIGHT  DATA  FROM  NAGRA  TAPE  TO  DEC  SYSTEM-10 
DISK  FILES 

Flight  data  are  recorded  by  the  data  acquisition  system  [4] 
on  magnetic  tape  using  a  NAGRA  tape  recorder.  The  data  are  in  serial 
form  on  two  parallel  tracks,  as  shown  in  Fig.  2.  The  system  allows, 
without  sub-multiplexing,  20  channels  of  data  to  be  recorded  in  12-bit 
form  at  a  sampling  frequency  of  60  Hz  for  each  channel.  However,  by 
sub-multiplexing  each  of  three  channels  to  record  two  separate 
quantities  at  30  Hz,  and  each  of  another  three  to  record  four  separate 
quantities  at  15  Hz,  32  data  source  outputs  are  actually  recorded. 
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In  order  to  process  the  data  on  the  DEC  system-10  computer 
at  ARL,  two  transcription  processes  are  performed.  Firstly,  using  a 
specially  built  transcriber  unit  [5],  the  data  are  converted  from 
twin-track  to  single-track  forir,  as  shovn  in  Fig.  3,  and  stored  on 
computer  compatible  7-track  magnetic  tape.  Secondly,  by  running  the 
program  EFREE,  each  flight  record  on  the  7-track  magnetic  tape  is 
converted  to  a  disk  file.  The  program  provides  an  automatic  file 
labelling  system  whereby  the  user  specifies  the  first  one  to  three 
characters  of  a  five  character  name,  with  the  remaining  characters 
forming  a  numerical  sequence  starting  at  one  (e.g.  if  the  first  two 
characters  are  specified  as  4T,  the  sequence  of  filenames  would  be 
4T001,  4T002,  4T003,...).  On  subsequently  running  EFREE  with  the 
same  initial  character  string,  existing  files  are  not  overwritten  and 
the  numerical  sequence  is  continued.  Where  checksum  or  parity  errors 
are  found,  a  file  is  still  created,  but  with  nil  blocks.  The  following 
example  +  shows  the  use  of  EFREE  when  using  a  disk  pack  to  store  the 
disk  files.  In  the  example,  before  running  EFREE,  the  7-track  magnetic 
tape  and  disk  pack  SKI  are  first  mounted  and  then  assigned  according 
to  the  requirements  specified  by  EFREE.  To  enable  the  files  to  be 
written  directly  onto  the  disk  pack,  the  users  disk  area  search  list 
must  first  be  modified  by  running  the  system  program  SETSRC.  Following 
conversion  of  the  last  complete  flight  record  on  the  7-track  tape, 
in  attempting  to  read  and  transcribe  another  record,  EFREE  usually 
terminates  because  of  an  error  condition. 


.PLE  CHECK  MX  MTA/TTR24 
MESSAGE  TRANSMITTED  at  11:41 
.PLE  CHECK  MOUNT  DSKPACK  SKI 
MESSAGE  TRANSMITTED  at  11:42 

• 

; ;OPR:  -  ON  MTA3 

• 

; ;OPR:  -  AND  SKI  IS  MOUNTED 

.AS  MTA3  MTI  Assignment  required  by  EFREE 

MTA003  assigned 


f  In  coaputer  program  examples,  terminal  responses  typed  by  the  user 
are  underlined.  Comments  are  shown  in  a  different  type  face  along¬ 
side  the  appropriate  terminal  message. 
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.REW  MT1  Rewind  magnetic  tape 

.AS  DSK  MXO  Assignment  required  by  EFREE 

DSK  assigned 


.R  SET SRC 


*T  Type  existing  search  list 

DSKC: ,  DSKB: ,  FENCE 

*C  SKI , DSKC  Add  SKI,  berore  DSKC 

*T  Type  new  search  list 

SKI:,  DSKC:,  FENCE,  DSKB: 

.RU  EFREE 

ASSIGN  INPUT  MAGTAPE  AS  "MTI"  AND  OUTPUT  MAGTAPE  OR 
DISK  AS  "KTO”. 

OK  TO  PROCEED?  TYPE  "C"  TO  CONTINUE,  "S"  TO  STOP  :  C_ 

TYPE  IN  UP  TO  3  CHAR  FOR  FILENAME  :  24  Files  24001  to  24013  written 

onto  disk  pack  SKI 

? 

?Hung  device  MTA003;  UUO  at  user  PC  006741  EFREE  terminated  by  error 

condition 


3 .  DESCRIPTION  OF  DATA  REDUCTION  PROCEDURE 

The  complete  data  reduction  procedure  may  be  divided  into  a 
number  of  separate  stages.  These  are  input  and  basic  processing, 
correction  for  drop-outs,  correction  for  instrument  and  analogue  filter 
time  delays,  conversion  to  engineering  units,  filtering,  smoothing, 
calculation  of  additional  quantities,  and  output.  Apart  from  system 
programs  (e.g.  PRINT,  PLOTQ) ,  the  two  computer  programs  REFINE  and 
TRANS  should  be  used  each  time  data  are  processed.  REFINE  reads  in 
and  processes  the  flight  data,  storing  the  reduced  data  in  binary 
form  suitable  for  input  to  the  general  purpose  output  program  TRANS 
[7,  8],  which  produces  tabular  and  graphical  output.  Another  program, 
called  LABCAL,  is  used  to  create,  or  modify  if  already  existing,  a  file 
containing  information  (e.g.  calibrations) ,  required  by  REFINE,  that 
specifies  how  each  channel  is  to  be  processed.  The  file  is  referred 
to  as  a  calibration  file.  For  checking  purposes,  a  further  program. 


called  DDMP12,  may  be  used  to  output  raw  flight  data  in  both  octal 
and  decimal  form.  In  this  section,  only  the  processing  stages 
performed  by  REFINE  are  outlined.  Examples  on  the  running  of  REFINE, 
as  well  as  the  other  three  programs,  form  the  contents  of  Section  4 
below. 


3.1  Input  and  Basic  Processing 

Data  are  read  in  the  form  of  36-bit  words,  each  word  combining 
three  successive  12-bit  words  recorded  on  NAGRA  tape  (see  Fig.  3) . 

A  complete  block  of  64  36-bit  words,  representing  data  for  eight 
complete  sampling  time  intervals,  are  read  into  a  buffer  at  a  time. 

If  a  read  error  occurs,  the  complete  block  is  replaced  by  the  block 
previously  read.  The  original  12-bit  words  are  then  extracted  from 
the  36-bit  words. 

The  flight  data  record  is  assumed  to  start  when  the  digital 
timing  word  (i.e.  clock  time)  of  consecutive  frames  is  sequential. 

This  condition  may  be  satisfied  by  any  of  the  four  possible  combinations 
of  clock  time  for  the  two  tracks  over  the  two  consecutive  time  frames 
(e.g.  time  given  by  Track  1  at  first  frame  and  time  given  by  Track  2 
at  second  frame) .  Previously  recorded  data  are  assumed  spurious  and 
therefore  ignored.  Following  the  start,  clock  time  is  defined  by 
incrementing  the  value  at  the  previous  frame.  If  the  actual  values 
read  for  both  tracks  differ  from  their  set  value,  a  time  error  is 
recorded . 


Signals  are  recorded  from  the  channels  in  a  regular  sampling 
order  at  equi-spaced  time  increments  throughout  each  complete  sampling 
time  interval  of  l/60th  second.  However,  it  is  assumed  that  all 
signals  recorded  within  this  sampling  interval  are  at  the  clock  time 
recorded  at  the  beginning  of  the  interval.  Hence,  a  maximum  relative 
time  delay  error  of  close  to  l/60th  second  is  tolerated. 

Values  for  channels  not  recorded  at  a  specific  time,  because 
of  sub-multiplexing,  are  obtained  by  linear  interpolation,  which  sllows 
a  common  sampling  frequency  of  60  Hz  to  be  assumed.  To  enable 
interpolation  for  channels  that  are  sub-multiplexed  at  the  higher  rate 
of  4  to  1,  values  at  up  to  three  time  increments  before  and  after  the 
current  time  increment  must  be  known.  Hence,  the  flight  record  is 
reduced  by  three  frames  at  both  beginning  and  end  (i.e.  zero  time  is 
at  the  fourth  frame) . 

The  actuax  clock  time  is  replaced  by  a  time  count  which  is 
set  to  zero  at  the  beginning  (i.e.  fourth  frame  provided  no  starting 
delay  is  specified  -  see  Section  3.8)  of  the  flight  record  and  is 
incremented  by  unity  at  each  sampling  interval.  However,  to  enable 
certain  checks  to  be  made  (see  Section  4.2),  the  five  least  significant 


I 
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digits  (in  octal)  of  the  clock  time  are  represented  directly  as  channel 
33  without  application  of  a  calibration  factor  and  offset. 

3.2  Correction  for  Drop-outs 

Having  first  obtained  a  flight  record  free  of  obvious  reading 
errors,  the  next  stage  is  to  correct  for  less  obvious  ones  caused  by 
misplaced  bits  in  the  12-bit  binary  words.  These  are  referred  to  as 
'drop-outs'  and  can  only  be  detected  when  one  of  the  most  significant 
bits  is  misplaced. 

For  a  sequence  of  signals  jcQ,  x^,  x2,  which  includes  inter¬ 
polated  values  for  sub-multiplexed  channels,  a  drop-out  is  considered 
to  have  occurred  when 


+  x. 


-  x 


where  kQ  and  kj  are  constants  determined  appropriately  by  inspection 
of  the  flight  data.  For  a  12-bit  word,  giving  an  integer  value,  or 
count,  in  the  range  0  to  4095,  values  found  generally  appropriate  were 
ko  =  100  and  kj_  ■  5.  When  a  drop-out  occurs,  the  current  value  is 
reset  to  its  previous  value,  i.e.  X2  =  Xj_.  The  constants  kQ  and  k^ 
are  stored  for  convenience  in  the  calibration  file  in  the  location 
labelled  CLOCK  TIME,  i.e.  channel  33.  It  should  be  noted  that  these 
calibrations  are  not  used  in  any  way  for  the  clock  time. 


3.3  Correction  for  Instrument  and  Analogue  Filter  Time  Delays 

Time  delays  requiring  correction  are  inherent  in  some 
instrument  outputs  and  are  also  introduced  by  the  use  of  analogue 
filters  in  the  data  acquisition  system.  Consistent  with  inaccuracies 
already  accepted  by  assuming  that  all  signals  are  at  the  recorded  clock 
time,  these  corrections  are  expressed  in  terms  of  complete  sampling 
time  intervals. 


At  present,  the  only  instruments  where  significant  inherent 
time  delays  are  apparent  in  the  signal  output  are  those  recording 
ground  speed  (i.e.  longitudinal  and  lateral  doppler  velocity).  As 
part  of  the  post-flight  kinematic  consistency  checking  process  [6J , 
comparison  of  doppler  velocity  measurements  with  those  obtained  on 
integrating  linear  accelerometer  measurements  indicate  a  time  delay 
correction  of  33  sampling  intervals  (i.e.  0.55s). 

Sixth-oraer  analogue  Butterworth  lowpass  filters  were  used 
on  most  channels  in  the  data  acquisition  system,  with  cut-off  frequencies 
of  3, 6, and  12  Ha,  generally  corresponding  to  the  sampling  frequencies 
of  15,30,  and  60  Hz.  A  description  of  the  quantities  measured,  together 
with  their  sampling  frequency  and  analogue  filter  cut-off  frequency. 
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are  given  in  Table  1.  For  the  purpose  of  obtaining  data  suitable  for 
validation  of  the  Sea  King  mathematical  model  [2,3],  accurate 
correction  of  phase  shift  is  considered  necessary  only  at  lower 
frequencies  in  the  passband  of  each  filter.  In  this  region,  the  filters 
used  have  a  reasonably  linear  phase  response  (Ref.  9,  p.112)  so  that 
a  constant  time  delay  may  be  assumed.  For  the  filters  with  cut-off 
frequencies  of  3,6,  and  12  if z,  the  corresponding  time  delays  are  12, 

6,  and  3  sampling  intervals. 

3.4  Conversion  to  Engineer inc  Units 

To  convert  each  recorded  output  signal  from  an  integer  value 
x  in  the  ranae  0  to  4095  to  a  corresponding  value  y  in  engineering 
units,  the  signal  is  first  multiplied  by  a  calibration  factor  c  and 
then  added  to  an  offset  value  d,  i.e. 

y  =  cx  +  d 

These  calibration  constants  were  obtained  in  the  first  place  from 
manufacturer' s  specifications  and  special  calibration  tests.  However, 
in  the  light  of  post-flight  analysis  (e.g.  kinematic  consistency 
checking) ,  some  adjustments  were  nade  (6) . 

3.5  Digital  Filtering 

Digital  filtering  procedures  have  been  developed  [10]  to 
remove  noise  arising  from  rotor-induced  vibration,  such  noise  being 
excessive  in  some  channels  despite  analogue  pre- filtering .  Digital 
Butterworth  lowpass  filters  are  used,  optionally,  with  characteristics 
selected  to  suit  the  particular  variable  to  be  filtered.  Phase  shift 
is  corrected  on  the  assumption  tint  it  is  proportional  to  frequency, 
i.e.  there  is  a  constant  time  shift.  Two  standard  fifth-order  filters 
have  been  found  adequate  for  all  channels.  One  has  a  high  attenuation 
of  50  dB  at  the  rotor  frequency  of  3.5  Hz,  giving  a  low  cut-off 
frequency  of  1.12  Hz,  and  removes  noise  at  and  above  the  rotor  frequency. 
It  is  suitable  for  variables  where  the  energy  is  confined  to  relatively 
low  frequencies.  The  ether  filter,  which  has  a  higher  cut-off  frequency 
of  4  Hz,  is  used  alternatively  where  the  high  frequency  content  is 
significant. 

3.6  Smoothing 

Another  means  of  removing  noise  [11]  has  been  incorporated  in 
the  data  reduction  procedure.  In  the  method,  referred  to  as  smoothing, 
each  observation  is  replaced  by  the  value  lying  on  a  polynomial  curve 
fitted,  using  the  least  squares  criterion,  through  a  band  of  observations 
centred  on  the  current  observation.  The  implementation  described  in 
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Ref.  1  l’t  allows  the  tire  consuming  operation  of  matrix  inversion  to 
bo  replaced  by  the  use  of  recursive  expressions.  However,  to  avoid 
a  significant  accumulation  of  truncation  errors  with  this  method,  the 
least  squares  parameters  and  covariance  matrix  are  recalculated  using 
matrix  inversion  after  every  150  observations.  As  with  filtering, 
smoothing  for  each  channel  is  optional,  with  appropriate  characteristics 
selected.  Though  intended  as  an  alternative  to  filtering,  the  data 
may  be  smoothed  following  filtering.  The  main  advantage  smoothing 
has  over  filtering  is  that  no  phase  shift  is  introduced.  However,  for 
the  purpose  of  removing  noise  at  known  frequencies,  filtering  has  been 
found  to  be  more  effective  and  less  time  consuming. 

3.7  Calculation  of  Additional  Quantities 


There  are  a  number  of  additional  quantities  that  may  be  derived 
from  the  quantities  actually  measured  in  the  data  acquisition  system. 

Some  of  the«e  quantities  are  required  for  validation  of  the  Sea  King 
mathematical  model,  while  others  are  necessary  for  confirmation  of 
the  consistency  of  the  data  itself  in  the  measurement  of  kinematic 
quantities.  The  additional  quantities  are  therefore  divided  into  four 
groups,  which  include  (1)  quantities  obtained  using  pressure  measurements 
from  a  boom-mounted  pitot-static  probe,  (2)  blade  anqles,  (3)  Euler 
angles,  and  (4)  quantities  for  kinematic  consistency  checking.  Table  2 
provides  a  list  of  each  quantity  belonging  to  the  above  four  groups. 
Documentation  on  the  calculation  of  the  quantities  is  given  in  Ref.  1 
for  the  first  two  groups  and  in  Ref.  6  for  the  last  two  groups. 

Table  3  summarises  the  use  made  in  the  calibration  file  of  the 
calibration  factors  and  offsets  for  these  additional  quantities 
together  with  the  treasured  quantities  and  clock  time. 

3 . G  Output 

The  output  of  all  channels  for  a  complete  flight  record  is 
often  not  required  for  reasons  of  economy  of  computing  time,  disk 
storage  space,  and  hard  copy  output.  The  program  REFINE  therefore 
allows  the  user  to  restrict  output  to  specific  channels,  either 
individually  or  as  groups,  and  to  a  specific  time  range  and  interval. 
Though  presenting  proportional  savings  in  disk  storage  space  and  hard 
copy  output,  these  restrictions  do  not  necessarily  allow  corresponding 
savings  in  computing  time.  The  reasons  for  this  are  as  follows. 

Firstly,  even  though  some  channels  may  not  be  required  for  output,  they 
may  still  be  required  to  be  calculated  in  order  to  derive  additional 
quantities  dependent  on  them.  Secondly,  for  an  output  interval  areater 
than  the  sampling  trequency,  it  is  still  necessary  to  calculate  quantities 
(a)  at  each  sampling  time  increment  when  filtering,  smoothing,  or 
correcting  for  time  delays,  and  (b)  at  up  to  three  time  increments  before 
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and  after  the  current  tire  increment  to  enable  interpolation  for 
channels  sub-multiplexed  4  to  1 .  It  should  also  be  noted  that  every 
frame  of  data  "ust  be  read  ir.  up  to  the  starting  tin j,  as  v/ell  as 
beyond  the  upper  tine  lir-.it  to  allow  for  the  maximum  accumulated  tine 
delay. 


Having  obtained  the  required  output  in  binary  form  from 
REFINE,  the  program  TRANS  is  then  used  to  convert  the  output  into 
tabular  and  graphical  torn.  because  TRAPS  is  also  used  in  the  same 
way  for  the  output  of  the  Sea  King  mathematical  model,  comparisons 
between  flight  data  and  model  prediction  on  the  same  graph (s)  are 
readily  made  using  the  'REPEAT1  command  facility.  In  doing  this, 
however,  assigned  block  numbers,  corresponding  to  those  used  in  the 
mathematical  model,  should  be  used  for  the  flight  data  instead  of 
channel  numbers. 

4.  EXAMPLES  OF  USE  OF  COMPUTER  PROGRAMS 


The  use  of  the  four  computer  programs  LAP.CAL,  DDI1P12,  REFINE, 
and  TRANS  is  now  illustrated  by  examples,  with  the  complete  processing 
procedure  divided  into  three  stages.  Firstly,  the  use  of  LABCAL  is 
demonstrated  in  creating  a  calibration  file  and  then  nodifyinq  it. 
Secondly,  initial  checks  on  the  correctness  of  the  data  reduction 
procedure  are  made  by  com^arino  uncaiibrated  data  in  integer  form 
given  directly  by  DDUP12  with  those  obtained  using  REFINE  and  TRAMS. 
Having  once  established  the  correctness  of  the  data  reduction  procedure, 
this  stage  is  generally  omitted  when  processing  further  flight  records. 
Thirdly,  subsequent  processing  is  performed,  again  using  REFINE  and 
TRANS,  in  which  various  forms  of  output  of  fully  processed  data  are 
produced . 

4.1  Creating  ojf  Modifying  Calibration  File 

The  program  LABCAL  has  boon  written  to  enable  the  calibration 
file  required  in  a  specific  format  by  REFINE  to  be  easily  created,  and 
then  modified.  In  making  modifications,  the  particular  variable,  or 
group  of  variables,  to  be  modified  is  specified  using  a  key  number. 

These  numbers  are  defined,  when  appropriate,  by  the  terminal  responses 
of  LABCAL.  The  following  example  shov;s  how  such  a  file  (Fig.  4) , 
restricted  here  to  four  channels  for  brevity,  may  first  be  created. 
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.RU  LABCAI. 

NUMBER  OF  CHANNELS  (excluding  TIME;  =  4 
CALIBRATION  FILENAME  "  CAL 4 
[Cre«cir.,T  new  calibration  file] 

[To  exit,  type  EXIT  frr  label;  unspecified  information 
then  stored  as  zeroes  or  blanks] 

[If  computer  crashes,  unspecified  information  stored  as 
zeroes  or  blanks;  re-run  and  modifv] 

[To  include  exceptions,  re— run  arm  modify] 

TITLE  (2  lines  of  60  cbrs) 

:TEST  CALIBRATION  Eli.E 

: 4  CHANNELS  ONLY 

TIKE  LABEL  (2  lines  of  !0  chrs) 

:  TIME 

: _ (S) 

TIME  CAL  FACTOR,  OFFSET  =  0.016667, 0 

[For  no  set  plot  limits,  answer  by  tynins  "c-r"} 

***  CHANNEL  i  *** 

LABEL  (2  line?  of  10  ch.-s) 

■  CLAN  1 

CAL  FACTOR,  OFFSET  -  1,0 
ASSIGNED  3LK  NUMBER  =  _i_ 

PLOT  LIMITS  fLOWER,  UPPER )  =  0, 10 

FILTER  FREQ  'Hz),  ATTENUATION'  fcS)  =  3.5,50 

NO.  POLES  «  2.  . 

NO .  PNTS  SMOOTHED ,  NO .  P ARAMS  -  __ 

***  CHANNEL  2  *** 

LABEL  (2  lines  of  10  chrs) 

;  CHAN  2 

CAL  FACTOR,  OFFSET  =  0. 12,-0. 4 
ASSIGNED  BLK  NUMBER  =  5 
PLOT  LIMITS  (LOWER,  UPPER )  =  _ 

FILTER  FREQ  (Hz),  5TTr'ii'  >ja>  = 

NO.  POLES  =  _ 

NO.  PNTS  SMOOTHED,  NO.  PARAMS  =  0,2 


***  CHANNEL  3  *** 

LABEL  (2  lines  of  10  chrs) 
:EXIT 


STOP 

END  OF  EXECUTION 

CPU  TIME:  2.00  ELAPSED  TIME:  2:46.38 
EXIT 
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In  the  above  example,  calibration  information  is  specified 
only  for  the  first  two  channels .  By  typing  ‘EXIT’  for  the  first  line 
of  the  label  of  channel  2,  the  unspecified  information  for  channels 
3  and  4  is  stored  as  zeroes  or  blanks.  To  complete  the  calibration 
file  (Fig.  5) ,  it  may  be  modified  by  again  running  LABCAL  as  follows. 

.RU  LABCAL 


NUMBER  OF  CHANNELS  (excluding  TIME'  =■  4 
CALIBRATION  FILENAME  =  CAL 4 
[Calibration  file  already  exists] 

ARE  MODIFICATIONS  REQRD  :  Y 

[Typo  Kev  no.;  =0  (all  vars);  *  i  (label  -  Inc.  title); 

=  2  ( cal  factor,  offset);  «  3  (assigned  no.); 

*  4  (plot  limits):  =  5  (filter  characteristics); 

=  6  (smoothing  characteristics);  =  7  (exceptions); 
c  8  (channel  no.)] 

[Channel  no.  -1  denotes  TIME,  -2  denotes  TITLE] 

KEY  NUMBER  ’  0 

CHANNEL  NUMBER  =  3 
LABEL  (2  lines  of  10  chrs) 

:  CHAN  3 

CAL  FACTOR,  OFFSET  =  8. 3,0.001 
ASSIGNED  BLK  NUMBER  *  24 
PLOT  LIMITS  (LOWER,  UPPER)  =  _ 

FILTER  FREO  (Hz)  ,  ATTENUATION  (dB)  =  4.0, 3.0 
NO.  POLES  =  5 

NO,  PNTS  SMOOTHED,  EC.  F/.f.AMS  -  9,3 

CHANNEL  NUMBER  =  4 
LABEL  (2  lines  of  10  chrs) 

:  CHAN  4 

CAL  FACTOR,  OFFSET  *  0.00254,-0.003 
ASSIGNED  BLK  NUMBER  =  2 
PLOT  LIMITS  (LOWER,  UPPER)  =  -2.2 
FILTER  FREQ  (Hz),  ATTENUATION  (dB)  =*  _ 

NO.  POLES 

NO.  PNTS  SMOOTHED,  NO.  P ARAMS  =*  _ 

CHANNEL  NUMBER  ■=  _ 

ARE  MODIFICATIONS  REQRD  :  N 
STOP 

END  OF  EXECUTION 

CPU  TIME:  1.50  ELAPSED  TIME:  8.24 
EXIT 
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Because  all  the  variables  were  required  to  be  specified  for 
channels  3  and  4  in  che  above  example,  a  key  number  of  0  was  used. 
Generally,  when  further  modifications  are  required,  one  variable  needs 
to  be  changed  for  one  or  more  channels.  For  key  numbers  1,2, 3, 5,  and  6, 
the  individual  terminal  messages  typed  by  the  program  are  identical  to 
those  given  in  the  above  example.  For  other  key  numbers,  i.e.  4,7,  and 
8,  further  explanation  is  now  given. 

For  key  number  4,  as  well  as  being  able  to  specify  plot 
limits  directly  as  above,  they  may  be  read  from  a  specified  disk  file, 
in  the  format  of  a  scale  limits  file  for  TRANS,  i.e.  TRANS. SCA  or 
TRANS. SCL  file.  This  facility  allows  plot  limits  that  have  been 
modified,  ox  created,  on  running  TRANS,  to  overwrite  the  limits 
originally  specified  in  the  calibration  file. 

Exceptions  for  calibration  factors  and  offsets  appropriate 
to  a  particular  flight  can  only  be  specified  by  using  key  number  7. 

The  exceptions  are  listed  at  the  end  of  the  calibration  file. 

It  may  be  necessary  to  insert  or  delete  quantities  in  the 
calibration  file,  thus  requiring  an  incremental  change  in  all  channel 
numbers  following  the  change,  key  number  8  allows  a  range  of  channel 
numbers  to  be  incremented  by  a  constant  amount,  positive  or  negative. 

New  quantities  may  then  be  inserted  using  a  test  editing  program  such 
as  TECO.  It  should  be  noted  that  when  using  key  number  8,  the  channel 
numbers  in  the  exceptions  list  are  unaffected  and  will  need  to  be  edited 
separately. 

The  following  example  illustrates  the  use  of  key  numbers 
4,7,  and  8,  with  the  resulting  calibration  file  shown  in  Fig.  6. 


.TYPE  TRANS. SCA 

1  0. 0000E+00  2. OOOOE+Ol 

3  -5 . 0000F.+Q0  5. 0000E+Q0 

•RU  LABCAL 


NUMBER  OP  CHANNELS  (excluding  TIME)  **  4 
CALIBRATION  FILENAME  -  GALA 
[Calibration  file  already  exists] 

ARE  MODIFICATIONS  REORD  :  Y 

[Type  Key  no.;  =0  (all  vars) ;  **  I  (label  -  inc.  title); 

**  2  (cal  factor,  offset);  «  3  (assigned  no.); 

*  4  (plot  limits);  *  5  (filter  characteristics); 

“  6  (smoothing  characteristics);  «  7  (exceptions); 
*=  8  (channeL  no.)] 

[Channel  no.  -1  denotes  TIME,  -2  denotes  TITLE] 


KEY  NUMBER  =  4 

ARE  PLOT  LIMITS  TO  HE  REAP  FROM  DS  K  :  N 

CHANNEL  NUMBER  =  2 

PLOT  LIMITS  (LOWER,  ITFF.R  )  =  -50,0 

CHANNEL  NUMBER  -  . 

ARE  MODIFICATIONS  REQRD  :  Y 

(Type  Key  no.;  =  0  (all  vars) ;  =  1  (label  -  Inc.  title); 

=*  2  (cal  factor,  offset);  =  3  (assigned  no.); 

■  4  (plot  Units);  =  5  (filter  characteristics); 

*  6  (smoothing  characteristics);  =*  7  (exceptions) 

*  8  (channel  no.); 

[Channel  no.  -1  denotes  TIME,  -2  denotes  TITLE] 

KEY  NUMBER  -  4 

ARE  PLOT  LIMITS  TO  BE  READ  FROM  DSK  :  Y 

FILENAME  -  TRANS . SCA 

ARE  ASSIGNED  BLK  NUMBERS  USED  :  N 

ARE  MODIFICATIONS  REQRD  :  Y 

[Type  Key  no.:  =  0  (all  vars);  =  1  (label  -  inc.  title); 

=*  2  (cal  factor,  offset);  =  3  (assigned  no.); 

=  4  (plot  limits);  -  5  (filter  characteristics); 

=  6  (smoothing  characteristics);  =  7  (exceptions) 
•»  8  (channel  no.)] 

[Channel  no.  -1  denotes  TIME,  -2  denotes  TITLF. ! 

KEY  NUMBER  “  7 

CHANNEL  NUMBER,  FLIGHT  NUMBER  «  4^2 
CAL  FACTOR,  OFFSET  -  0.0254,-0.015 

CHANNEL  NUMBER,  FLIGHT  NUMBER 

ARE  MODIFICATIONS  P.EQRD  ;  Y 

[Type  Key  no.;  =0  (all  vars);  *  1  (label  -  inc.  title); 

“  2  (cal  factor,  offset);  *  3  (assigned  no.); 

=  4  (plot  limits);  =  5  (filter  characteristics); 

=  6  (smoothing  characteristics);  **  7  (exceptions) 

*  8  (charnel  no-)] 

[Channel  no.  -1  denotes  TIME,  -2  denotes  TITLE] 

KEY  NUMBER  -  8 

CHANNEL  NUMBER  LIMITS;  LOWER,  UPPER  =  3, 4 
ADJUSTMENT  FOR  EACH  CHANNEL  -  2 

ARE  MODIFICATIONS  REQRD  :  N 

STOP 

END  OF  EXECUTION 

CPU  TIME:  2.70  ELAPSED  TIME:  15.78 
EXIT 
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4.2  Initial  Chock  on  Uncalibrated  Data 


In  order  to  be  able  to  check  each  stage  of  the  data  reduction 
procedure  for  correctness,  a  iarqe  degree  of  flexibility  has  been 
incorporated  into  REFINE.  As  a  starting  pornt  in  the  checking  process, 
the  program  DDMP12  is  first  used  to  obtain  a  tabulation  (or  dump)  of 
unprocessed  flight  data  in  both  octal  and  equivalent  12-bit  decimal 
words.  It  should  be  noted  that  DDMP12  operates  in  terms  of  128 
36-bit  word  blocks,  which  are  twice  the  size  of  the  previously  defined 
flight  data  blocks.  Hence  each  block  dumped  will  represent  data  for 
16  time  increments.  For  flight  record  55025,  a  dump  of  the  first 
block,  given  in  Fig.  7,  is  obtained  as  follows. 


.RU  DDMP12 


Input  File  ?  :  15025 

Block  numbers  you  want  dumped  ?  (Specify  range,  e.g.  4-7)  :  1-1 
Output  File  ?  :  15025. DMP 

END  OF  EXEC OT ION 

CPU  TIME:  0.79  ELAPSED  TIME:  44.98 
EXIT 


Each  line  of  the  dump  represents  one  frame  or  time  increment. 
Individula  words  may  be  identified  using  the  single-track  form 
description  in  Fig.  3  md  the  sub-multiplexing  key  in  Fig.  2.  The 
program  REFINE  is  then  used  with  appropriate  options  selected  for 
minimal  processing  to  enable  a  comparison  of  results,  given  on  running 
TRANS,  with  those  already  given  by  the  dump.  The  new  results  will  of 
course  be  formatted  differently,  with  interpolated  values  shown  for 
channels  not  recorded  -*t  c.  specific  time  because  of  sub-multiplexing. 

The  calibration  file  CAL  used  when  running  REFINE  for  all  flight  data 
is  shown  in  Fig.  8.  However,  for  minimal  processing,  only  the  labelling 
information  is  actually  required  from  this  file. 

The  following  exarrple  shows  the  use  of  REFINE  and  TRANS  in 
producing  minimally  processed  output,  given  in  Fig.  9,  for  selected 
channels  at  the  beginning  of  flight  record  15025.  Clock  time,  represented 
by  channel  33  as  the  five  least  significant  digits  in  octal,  allows  ready 
comparison  between  Fig.  7  and  9.  Because  the  calibration  factors  and 
offsets  are  ignored,  the  units  given  in  the  channel  labels  are  inappropriate. 
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In  the  example,  it  is  assumed  that  the  flight  data  file  needs 
transferring  from  a  specific  area  on  the  disk  pack  SKI  to  the  user's 
own  disk  area . 


.PLE  CHECK  MOUNT  DSKPACK  SKI 
MESSAGE  TRANSMITTED  at  16:37 


; ;OPR:  -  OK  MOUNTED 
•  R  SETSRC 


*1 

DSKC:,  DSKB:  ,  FENCE 
*A  SKI: /NOWRITE 


or 


r£ 


COPY  /X -SKI: [1031 


1063] 15025 


*1 

DSKC:,  DSKB:,  SKI: /NOWRITE,  FENCE 
*"0 


.COPY  /X-[1031. 10631 15025 


.RU  REFINE 


INPUT  DATA  FILENAME  =  15025 

".DAT"  OUTPUT  FILENAME  (w/o  ext)  15025.DAT  by  default 

TITLE  (2  lines  of  60 

.  15025  by  default 

:RAW  DATA  IN  COUNTS  -  NO  DELAYS 
ARE  ASSIGNED  BLK  NUMBERS  REQRD  :  N_ 

CALIBRATION  FILENAME  =  CAL 
FLIGHT  NUMBER  -  3 

DOES  CHANNEL  18  MEASURE  TORQUE  :  N_ 

OUTPUT  INTERVAL  (in  60'ths  of  sec;  e.R.  12  for  0.2  sec)  =*  l 
STARTING  TIME  DELAY,  NO.  BLKS  FIRST  IGNORED  -  0,0 
TIME  LIMIT  -  0e3 
IS  FILTERING  REORD  :  N 
IS  SMOOTHING  REQRD  :  N 

ARE  INSTRUMENT  &  ANALOGUE  FILTER  DELAY  ADJUSTMENTS  REQRD  :  N 

ARE  SCALES  AND  OFFSETS  REQRD  :  N 

ARE  PLOT  LIMITS  REQRD  :  N 

ARE  DROP-OUTS  TO  BE  CORRECTED  :  N 
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ARE  ALL  CHANNELS  REQRD  :  N 
CHANNEL  NUMBERS  RF.QRP  FOR  FOLLOWING  CROUPS 
[Set  first  value  -ve  if  excluding] 

[Set  first  value  to  -100  if  non-2  excluded} 


Instumentation  data 
Boom  calculations 
Blade  angles 
Euler  angles 
Kinematic  consistency 
DATE  «  2-Apr-82 

TIME  =  13:  19 
No.  time  corrections  =* 
No.  blocks  replaced  = 

No.  drop-out  corrections 
STOP  000001 


(  1  to 
( 34  to 
(41  to 
( 4  6  to 
(48  to 


0 


33) 

40) 

45) 

47) 

70) 


13,17,18,21,23,26.28.33 


Error  corrections  in  15025. ERR 


END  OF  EXECUTION 

CPU  TIME:  3.98  ELAPSED  TIME:  23.24 
EXIT 


.RU  PUB: TRANS 


[TRANS  version  date  12-NOV-801 

I/P  FILENAME  -  15025 

15025 

RAW  DATA  IN  COUNTS  -  NO  DELAYS 

I/P  FILE  RECORDED  ON  2-Apr-82  AT  13:19 

INTEGN  INT  -  O.OOOOE+OO;  RUN  CPU  TIME  =  0.00  SEC. 

TIME  FROM  0. 0000F.-K)0  TO  3.0000E-01  IN  STEPS  OF  1.6667E-02 
*PRC 

PRINTING  IN  COLUMNS  : 

BLKS 

A. 

IS  0/P  TO  TTY  RF.QRD  :  N 
*GOE 

**  RUNNING  ** 

*EXI 

END  OF  EXECUTION 

CPU  TIME:  1.12  ELAPSED  TIME:  2.70 
EXIT 


i. 
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Having  verified  the  correctness  of  the  minimally  processed 
data  in  Fig.  9,  the  various  stages  in  the  data  reduction  procedure 
may  be  verified  one  at  a  time.  Though  such  checks  have  been  performed, 
only  the  results  of  the  next  stage  of  the  process  are  presented. 

Fig.  10  shows  the  data  as  in  Fig.  9,  except  that  instrument  and 
analogue  filter  time  delays  are  included.  The  clock  time  is  now 
correctly  referenced  only  to  the  channel (s)  with  the  maximum  time 
delay.  In  this  case,  channel  21  has  a  total  time  delay  of  39  sampling 
increments  (0.65s) ,  which  means  adding  47  (in  octal)  to  the  clock  time 
in  Fig.  9  to  arrive  at  the  value  in  Fig.  10. 


4.3  Subsequent  Proces sing 


Assuming  that  all  stages  of  the  data  reduction  procedure  have 
been  checked  for  correctness,  REFINE  may  now  be  used  with  confidence  to 
obtain  fully  processed  data.  The  following  exanqple  illustrates  many  of 
the  capabilities  of  REFINE  used  together  with  TRANS  in  producing  tabular 
and  graphical  output.  Files  containing  filter  characteristics  and 
velocity  iteration  information  are  produced  by  REFINE  in  the  example. 

The  user  should  refer  to  Refs.  10  and  6  respectively  for  further  information 
on  the  quantities  printed.  An  additional  34  (in  octal)  should  now  be 
added  to  the  clock  time  in  Fig.  10  because  of  the  additional  time  delay 
of  28  sampling  increments  (0.467  s)  due  to  digital  filtering  (actual 
delay  calculated  is  0.459  s)  . 


INPUT  DATA  FILENAME  »  15025 
".DAT"  OUTPUT  FIt.ENAMF.  (w7  o  ext)  ” 
TITLE  (2  lines  of  60  chrs) 


ARE  ASSIGNED  BLK  NUMBERS  REQRD  :  N. 

CALIBRATION  FILENAME  -  CAL 
FLIGHT  NUMBER  -  3 

DOES  CHANNEL  18  MEASURE  TORQUE  :  N 

OUTPUT  INTERVAL  (in  CO'^ho  o L  sec;  e.g.  12  for  0.2  sec)  ■  6_ 

STARTING  TIME  DELAY,  NO.  BLKS  FIRST  IGNORED  -  0,0 

TIME  LIMIT  *  8.  Zero  or  blank  gives  complete  record 

IS  FILTERING  REQRD  :  Y 

IS  OUTPUT  OF  FILTER  CHARACTERISTICS  REQRD  :  Y 
ARE  DIGITAL  FILTER  DELAY  ADJUSTMENTS  REQRD  ;  Y 
IS  SMOOTHING  REQRD  :  N 

ARE  INSTRUMENT  &  ANALOGUE  FILTER  DELAY  ADJUSTMENTS  REQRD  :  Y 

ARE  SCALES  AND  OFFSETS  REQRD  :  Y 

ARE  PLOT  LIMITS  REQRD  :  Y 

ARE  DROP-OUTS  TO  BE  CORRECTED  :  N 

ARE  ALL  CHANNELS  REQRD  :  Y 

DATE  “  5-Apr-82 

TIME  -  09:39 

IS  VELOCITY  ITERATION  INFORMATION  REQRD  :  Y 

ARE  ANY  FAULTY  DATA  TO  BE  REPLACED  BY  ALTERNATIVELY  DERIVED  DATA  :  U 
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No.  time  corrections  =  0 

No.  blocks  replaced  --  0 

No.  drop-out  corrections  =  0 

STOP  000001 

END  OF  EXECUTION 

CPU  TIME:  1:33.78  ELAPSED  TIME:  4:10.26 

EXIT 

.RU  PUB: TRANS 

£ 


[TRANS  version  date  12-NOV-80] 

I/P  FILENAME  =  15025 
15025 

I/P  FILE  RECORDED  ON  5-Apr-82  AT  09:39 

INTEGN  INT  -  0.0000E+00;  RUN  CPU  TIME  =  0.00  SEC. 

TIME  FROM  0.  OOOOE+OO  TO  7.4000F.400  IN  STEPS  OF  1.0000E-01 
*SCA 

BLK  NO.  -1  DENOTES  INDEP  VARIABLE 

ARE  PLOT  SCALE  LIMITS  TO  BE  READ  FROM  DSK  :  Y 

IS  TTY  LISTING  OF  LIMITS  REQRD  :  N 

ARE  MODIFICATIONS  REQRD  :  Y 

BLK,  LOWER,  UPPER 

15,-20,20 


4*7  • 

37.1900.2100 
38,60,100 
39,0,20 
50,-10.10 

54, -10, 10 

55, -10, 10 
57,-10,10 

56, -10, 10 

60.40.80 

63.40.80 

64.1900.2100 
67,-1 0,10 

*LI5T 

BLOCK  LIST  NUMBER  -  1 


LIST  OF  BLOCK  NUMBERS  -  1,2.3.16,5.6,7.27 
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BLOCK  LIST  NUMBER  -  2 

LIST  OF  BLOCK  NUMBERS  »  8,11,17,15,10,18,29,30 
BLOCK  LIST  NUMBER  =  3 

LIST  OF  BLOCK  NUMBERS  -  12,13, 14,41,42,45,44 
BLOCK  LIST  NUMBER  -  4 

LIST  OF  BLOCK  NUMBERS  =  21,22,37,38,39,40,4,9 
BLOCK  LIST  NUMBER  =  _5 

LIST  OF  BLOCK  NUMBERS  -  54,55,56,57,58,59,62,63 
BLOCK  LIST  NUMBER  «  _ 

*TIM 

TIME  PARAMS;  LOWER,  UPPER,  INTERVAL  -  0,8,0. 5 
*PRC 

PRINTING  IN  COLUMNS  : 

BLKS 

13,17,18,21,23,26,28,33 

IS  O/P  TO  TTY  REQRD  :  N 
*SCA 

BLK  NO.  -1  DENOTES  INDEP  VARIABLE 

ARE  PLOT  SCALE  LIMITS  TO  BE  R^AD  FROM  DSK  :  Y 

IS  TTY  LISTING  OF  LIMITS  REQRD  ;  N 

ARE  MODIFICATIONS  REQRD  :  N 

*SPA 

IS  SPACING  BETWEEN  PLOTS  REQRD  :  N 

*PLS  Output  to  15025. PLT  (see  Fig.  14) 

STRIP  PLOTS  : 

BLKS 

U 

TO  SPECIFY  NO.  OF  X  UNITS /INCH,  TYPE  0  FOR  X 
LENGTH  OF  AXES  IN  INCHES;  X,  Y  -  0, 1 


Repeated  for 
lists  L2  - 


k 


block 

L5 


1 
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MIN  X,  NO.  OF  X  UNITS /INCH  =  0^ 

ARE  SYMBOLS  REQRD  FOR  PLOTS  :  N 
LINE  KEY  (0  GIVES  PFFA'F  T  '  0 

*GOE 

**  RUNNING  ** 

I 

I 

*EXI 

END  OF  EXECUTION 

CPU  TIME:  22.8.*  ELAPSED  II ME:  4;1.94 
EXIT 


For  kinematic  consistency  checking,  with  alternatively 
derived  quantities  shown  on  the  same  graphs,  TRANS  may  be  run  again 
as  follows. 


.RU  PUB: TRANS 


[TRANS  version  date  12-NOV-80] 

I/P  '"'LFMAME  -  1502, 

15025 

I/P  FILE  RECORDED  ON  5-Apr-82  AT  10:04 

INTEGN  INT  *  0.0000E-K50;  RUN  CPU  TIME  =*  0.00  SEC. 

TIME  FROM  0. 0000E+00  TO  7.4000E-+00  IN  STEPS  OF  1.0000E-01 
*SCA 

BLK  NO.  -1  DENOTES  INDEP  VARIABLE 


ARE  PLOT  SCALE  LIMITS  TO  BE  READ  FROM  DSK  :  Y 
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IS  TTY  LISTINC  OF  LIMITS  REQRD  :  N 
ARK  MODIFICATIONS  REQRD  :  N 

*SPA 

IS  SPACING  BETWEEN  PLOTS  REQRD  :  N 

*pL0  See  Fig.  15 

(PLS/O  Output,  for  this  run,  going  to  DSK:  15025.  PLO] 


OVERLAY  PLOTS  : 

BLKS 

46,68 

BLKS  TO  DETERMINE  Y  SCALE 
46 

TO  SPECIFY  NO.  OF  X  UNITS/INCH,  TYPE  0  FOR  X 

LENGTH  OF  AXES  IN  INCHES ;  X,  Y  «=  0^2 

MIN  X,  NO.  OF  X  UNITS/INCH  =  <^2 

ARE  SYMBOLS  REQRD  FOR  PLOTS  :  N 

LINE  KEY  (0  GIVES  DEFAULT)  =  0 

*GOE 

**  RUNNING  ** 


Repeated  for  blocks 

15.69 

47.70 
11,51 

8,52 

17,53 

21,60 

22,61 

37,64 


*EXI 


END  OF  EXECUTION 

CPU  TIME:  30.64  ELAPSED  TIME:  6:51.10 
EXIT 


On  subsequently  re-running  REFINE  with  the  same  ".DAT" 
output  filename,  the  existing  filenames  15025. DAT  and  15025. ERR  are 
first  renamed  to  OLD. DAT  and  OLD. ERR  respectively. 
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5 .  CONCLUDING  REMARKS 

Data  reduction  procedures  used  in  obtaining  fully  processed 
data  from  raw  flight  data  for  trials  on  a  Sea  King  Mk.50  helicopter 
have  been  presented.  The  trials  were  conducted  for  the  purpose  of 
validating  the  Sea  King  mathematical  model.  The  procedures  may  be 

divided  into  a  number  of  stages,  which  allow  various  corrections  and  < 

calibrations  to  be  applied,  removal  of  noise  using  digital  filtering  ' 

or  'least  squares'  smoothing,  and  calculation  of  many  additional  > 

quantities.  Some  of  these  quantities  are  required  for  validation  of  ] 

the  mathematical  model,  while  others  are  necessary  for  confirmation 
of  the  consistency  of  the  data  itself  in  the  measurement  of  kinematic 

quantities.  The  use  of  the  various  computer  programs  developed  has  , 

been  illustrated  by  examples.  The  two  principle  programs,  REFINE  and 

TRANS,  should  be  used  each  time  data  are  processed.  REFINE  reads  in 

and  processes  the  flight  data,  storing  the  reduced  data  in  binary 

fom  suitable  for  input  to  the  general  purpose  output  program  TRANS, 

which  produces  tabular  and  graphical  output.  V.’ith  REFINE,  the 

facility  for  excluding  one  or  more  processing  stages  has  been  found 

very  useful  in  checking  the  correctness  of  each  of  these  stages. 

TRANS,  which  is  also  used  to  obtain  output  for  the  mathematical  model, 
has  been  found  sufficiently  versatile  in  providing  multi-curve  plots 
on  the  same  graph.  These  plots  ray  represent  more  than  one  variable 
in  a  single  flight  data  record,  as  required  when  checking  for  kinematic 
consistency,  or  the  same  variable  from  a  number  of  different  files, 
as  when  comparing  flight  trials  data  with  mathematical  model  results. 
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TABLE  1 

Quantities  measured  by  data  acquisition  system 


Channel  Quantity  measured  Symbol  Sampling  Analogue  filter 

number  freq. (Hz)  cut-off  freq. (Hz) 


2 

3 

4 

5 

6 

7 

e 

9 

10 

11 

12 

13 

14 

15 

16 
17 
18* 

19 

20 
21 
22 

23 

24 

25 

26 
27 


Cyclic  stick  position  -  pitch 
Cyclic  stick  position  -  roll 
Collective  stick  position 
Angle  of  attack 

Fore--aft  (pitch)  push-pull  rod 
position 

Lateral  (roll)  push-pull  rod 
position 

Collective  (altitude)  push-pull 
rod  position 


stk 

’’stk 


Cstk 

a 


d 

auxp 


d 

auxr 


d 

auxa 


Pitch  rate 
Sideslip  angle 
Roll  attitude 
Roll  rate 

Longitudinal  acceleration 
Lateral  acceleration 
Vertical  acceleration 
Pitch  attitude  (Euler  ancle) 
Yaw  pedal  position 
Yaw  rate 

Yaw  attitude  (heading) 

Lateral  cable  angle 
Longitudinal  cable  angle 
Doppler  lonqitudinal  velocity 
Doppler  lateral  velocity 
Doom  probe  dynamic  pressure 
Radio  altitude  (raw) 

Radio  altitude  ( smooth) 

Boom  probe  absolute  pressure 
Yaw  push-pull  rod  position 
Ambient  temperature 


g 

5 

P 

alonc 

alat 

a 

vert 

6 

d  . 
ped 

r 


0 

cH 

u' 

v' 


d 

auxy 

T  u 

amb 


•30 

60 

60 

60 

60 

GO 

60 

60 

60 

60 

60 

60 

60 

60 

30 

30 


30 


30 

30 

30 

15 

15 

15 

15 

15 

15 

15 


12 

12 

12 

12 

3* 

12 

12 

12 

12 

12 

12 

12 

12 

12 

6 

6 

6 


6 

6 

6 

6 

12 

3 

3 

3 

12  + 


28 


15 


i 


TABLE  1  (COtJT.) 


Channel 

number 

Quantity  measured 

Symbol 

Sampling 
freq. (Hz) 

Analogue  filter 
cut-off  freq. (Hz) 

29 

Torque  -  Engine  1 

Q1 

15 

- 

30 

Rotor  r . p . m . 

15 

- 

31 

Towed  probe  dynamic  pressure 

15 

3 

32 

Boom  probe  -  towed  probe 
differential  pressure 

Ap 

15 

3 

33 

Clock  time  -  in  octal  (5  least 
signif.  digits) 

tclock 

50 

- 

*  Switch  selectable  alternative  measurement  of  'Torque  -  Engine  2'  (Q^) 
t  Filters  inadvertently  interchanged. 


Notes:  1.  Push-pull  rod  positions  represent  auxiliary  servo 

displacements . 


2.  Acceleration  measurements  include  gravitational  effects. 

3.  Pitch  and  roll  attitude  are  relative  to  the  vertical 
(hence  roll  attitude  is  not  an  Euler  angle) . 

4.  Linear  accelerations  and  angular  velocities  are  in  body 
axes. 


5.  Doppler  velocities  are  in  moving  earth  axes. 

6.  Cable  angles  are  measured  with  respect  to  the  funnel 
axis  of  symmetry. 


I 
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TABLE  2 

Additional  quantities  calculated 


Channel 

number 

Quantity  calculated 

Symbol 

Group 

34 

Airspeed  from  boom  probe  at  standard  sea 
level  conditions,  corrected  for  position 
error 

V 

c 

1 

35 

Altimeter  setting  (12J 

oi:h 

1 

36 

Sea  level  temperature 

T 

o 

1 

37 

Altitude  from  boom  probe  static  pressure 

hb 

1 

38 

True  airspeed  at  aircraft  altitude 

VA 

1 

39 

Wind  velocitv 

vw 

1 

40 

Direction  from  which  wind  is  coming 

1 

41 

Longitudinal  cyclic  blade  pitch  angle 

Bls 

2 

42 

Lateral  cyclic  blade  pitch  angle 

Als 

2 

43 

Collective  blade  pitch  annle 

ec 

2 

44 

Tail  rotor  collective  blade  pitch  angle 

6T 

2 

45 

Collective  blade  pitch  anole  at  75%  rotor 

radius  position 

6C75 

2 

46 

Yaw  Euler  angle 

* 

3 

47 

Roll  Euler  angle 

$ 

3 

48 

Yaw  Euler  angle  derivative  -  by 
differentiation 

< ip 

D 

4 

49 

Fitch  Euler  anale  derivative  -  by 
differentiation 

e 

D 

4 

50 

Roll  Euler  angle  derivative  -  by 
dif  ferentiation 

4> 

D 

4 

51 

Roll  rate  -  by  differentiation 

PD 

4 

52 

Pitch  rate  -  by  differentiation 

qD 

4 

53 

Yaw  rate  -  by  differentiation 

rD 

4 

54 

Longitudinal  inertial  acceleration 

a 

X 

4 

55 

Lateral  inertial  acceleration 

a 

y 

4 

56 

Vertical  inertial  acceleration 

a 

z 

4 

57 

Longitudinal  velocity  derivative 

u 

4 

58 

Lateral  velocity  derivative 

V 

4 

59 

Vertical  velocity  derivative 

w 

4 

60 

Longitudinal  velocity  -  by  integration  (in 
moving  earth  axes) 

4 

TABLE  2  (CONT.) 


Channel 

number 

Quantity  Calculated 

Symbol 

Group 

61 

Lateral  velocity  -  by  integration  (in  moving 

earth  axes) 

v' 

I 

4 

62 

Vertical  velocity  -  by  integration  (in 

moving  earth  axes) 

Horizontal  velocity  magnitude  (u'j2+v'2)^ 

w 1 

I 

4 

63 

V7 

i 

4 

64 

Height  -  by  integration 

hi 

4 

65 

Yaw  Euler  angle  derivative 

4> 

4 

66 

Pitch  Euler  angle  derivative 

e 

4 

67 

Roll  Euler  angle  derivative 

i 

4 

68 

Yaw  Euler  angle  -  by  integration 

ip 

i 

4 

69 

Pitch  Euler  angle  -  by  integration 

e 

i 

4 

70 

Roll  Euler  angle  -  by  integration 

♦i 

4 

Notes:  1.  Linear  accelerations,  linear  velocity  derivatives,  and 

angular  velocities  are  in  body  axes,  unless  otherwise 
stated. 


2.  Euler  angles  are  measured  with  respect  to  earth  axes. 
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TABLE  3 

Use  of  calibration  factors  and  offsets 


Channel 
number (s) 

Use 

1-32 

Standard,  i.e.  y  =  cx  +  d 

33 

Drop-out  correction  constants  (see  Section  3.2)  with 

c  =  k  ,  d  =  k 
o  1 

34 

Airspeed  given  by  y  =  c  y(23)  +  d  where  dynamic 

pressure  q,  =  y(23) 
b 

35 

Altimeter  setting  given  by  y  =  c  =  QNH 

36 

Sea  level  temperature  given  by  y  =  c  =  Tq 

37 

d  =  offset  for  altitude  obtained  from  boom  probe 
static  pressure 

38  -  53 

No  use  made,  i.e.  c  =  d  »  0 

54  -  56 

d(54)  -  d(56)  give  offset  displacement  coordinates  of 
accelerometers  from  c.g.  position  [6] 

57  -  70 

Mo  use  made,  i.e.  c  *  d  *  0 

Mote:  x  a  output  signal  (integer  value  0  -  4095) 
y  =  output  signal  (engineering  units) 
c  «=  calibration  factor 
d  *  offset 

Muaiber  in  brackets  following  any  of  above  symbols  specifies 
particular  channel. 
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Axes  x  y_  z„ 
-  E  E  E 

x'  y 1  z' 


earth  (inertial)  axes 
moving  earth  axes 


x  y  z 


body  (helicopter  fuselage)  axes 


a)  Moving  earth  axes  relative  to  earth  b)  Body  axes  relative  to 
axes  (earth  axes  translated  such  moving  earth  axes 

that  origin  is  at  c.g.) 


Notes; 

1.  All  axes  systems  are  right-handed  orthogonal. 

2.  Horizontal  x'  axis  always  lies  in  fuselage  longitudinal  plane 
of  symmetry  (positive  forward) . 

3.  Moving  earth  axes  coincide  with  earth  axes  initial  time. 

4.  Orientation  of  body  axes  relative  to  earth  axes  are  defined  in 
conventional  way  for  aircraft,  i.e.  sucessive  rotations  in  yaw 

pitch  0,  and  roll  $.  Moving  earth  axes  represent  first 
intermediate  stage,  i.e.  rotation  in  yaw  ip. 


FIG.  1  Definition  of  axes  systems  and  their  relative  orientation. 


0 


Inter-frame  gap  - 
(1  x  12-bit  word) 


1  frame  =  12  x  12-bit  words 


NAGRA  tape 
(2  data  tracks  + 
voice  track) 


T  T 

LI  Ml 


Word  no. 


3  5  7  9  11  13 


3  4  5  6 


10  11 


Track  1 


T  T 
L2  M2 


4  6  8  10  12  14 


1812^^331, 


Track  2 


Notes: 

1.  Numbers  shown  in  each  word  box  are  transducer  channel  numbers 
(i.e.  after  sub-multiplexing) . 

2.  First  2  words  of  any  frame  on  each  track  form  a  digital  timing 
word  when  reversed  in  order  (T  is  'least'  signif.  word,  T^ 
is  'most'  signif.  word,  both  for  Track  1). 

3.  Last  3  words  of  any  frame  on  each  track  are  sub-multiplexed 
either  2:1  or  4:1  depending  on  2  least  significant  bits  of 
digital  timing  word  as  follows :- 


Timing 

bits 


Least  signif.  timing 

digit  (octal) 

0  or  4 

1  or  5 

2  or  6 

3  or  7 


Channels 

recorded 

Track 

1 

Track  2 

15 

19 

25 

17 

21 

29 

16 

20 

26 

18 

22 

30 

15 

19 

27 

17 

23 

31 

16 

20 

28 

18 

24 

32 

FIG.  2  Format  of  data  as  recorded  on  NAGRA  tape  recorder. 


Twin-track  form  (on  NAGRA  tape) 

1  frame  =  2  tracks  of  12xl2-bit  words 


ki 

T 

Ml 

i 

3 

5 

7 

• 

11 

13 

15 

17 

19 

TL2 

T 

M2 

2 

4 

_ 

6 

8 

10 

— 

E 

14 

16 

18 

20 

V 

Single-track  form  (on  7-track  magnetic  tape) 
1  frame  *  8x36-bit  words  =  24xl2-bit  words 


T  „ 

T  „ 

T 

R 

R 

L2 

Ml 

M2 

H 

■ 

10 


11 


12 

13 

a 

15 

16 

3 

18 


19 


20 


Note:  Numbers  shown  in  each  word  box  are  channel  numbers  before 
sub-multiplexing . 


FIG.  3  Conversion  of  data  from  twin-track  to  single-track  form. 
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FIG.  4  Newly  created  calibration  file  CAL4 
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FIG.  5  Initially  modified  calibration  file  CAM 
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FIG.  6  Subsequently  modified  calibration  file  CAL4 
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FIG.  7  File  15025. DMP  giving  dump  of  first  block  of  flight  record  15025 
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FIG.  8  Calibration  file  CAL  used  fox  all  flight  data 
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FIG.  9  File  15025. COr,  showing  minimally  processed  data  for  selected  channels  at  beginning  ot  flight  recce  d  1  MV 
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FIG.  13  File  1502 5. COIi  showing  fully  processed  data 
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FIG.  14a  Fully  processed  data  for  block  list  1 
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FIG.  14b  Fully  processed  data  for  block  list  2 
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FIG.  14e  Fully  processed  data  for  block  list  5 
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FIG.  15d  Kinematic  consistency  checking  for  longitudinal  and  lateral  velocity 


FIG.  I5e  Kinematic  consistency  checking  for  height 
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